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Abstract 

It has been some time since an empirical relation, which correlates DC with AC conductivity and 
contains a loosely defined coefficient thought to be of order one, was introduced by Barton, Nakajima 
and Namikawa. In this work, we derived this relation assuming that the conductive response consists of a 
superposition of DC conductivity and an AC conductivity term which materialized through a Havriliak- 
Negami dielectric function. The coefficient was found to depend on the Havriliak-Negami shape param- 
eters as well as on the ratio of two characteristic time scales of ions motion which are related to ionic 
polarization mechanism and the onset of AC conductivity. The results are discussed in relation to other 
relevant publications and they also applied to a polymeric material. Both, theoretical predictions and 
experimental evaluations of the BNN coefficient are in an excellent agreement, while this coefficient shows 
a gradual reduction as the temperature increases. 

Keywords: AC conductivity, DC conductivity, ionic conductors, ionic polarization mechanism, BNN 
relation 



1 Introduction 



Nearly four decades ago, an empirical relation was introduced by Barton, Nakajima and Namikawa which 
is known as BNN relation [1] [5J |3] . This expression correlates the electrical conductivity to the dielectric 
strength of the lower frequency polarization mechanism through, 

CTO = pSoAeUJmax , (1) 
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where p is a loosely defined parameter, expected to be of order 1, ctq is the apparent DC conductivity, 
Ae is the dielectric strength, uJmax is the angular frequency which corresponds to the maximum value of 
dielectric losses and Eq is the permittivity of vacuum. The loss peak which associated to Ae is, in general, 
characterized as broad. Ae may arise entirely from mobile charge effects and not involve bulk dielectric 
effects at all. When Eq. ([IJ is satisfied, both AC and DC conductivity may arise from the same type of 
charge transport mechanism [H [S] . 

The BNN relation has played an important role in the analysis and the treatment of frequency response 
dielectric data [SJ [TJ ISl HI UHl E] , assuming that p ^ \ with most of the corresponding works dealing with the 
scaling and universality issues of AC conductivity. In the literature a large number of different conductive 
disordered materials have been found to satisfy Eq. ([T]). The BNN relation is valid not only in amorphous 
solids, ionic glasses, single crystals and polymers but also in a variety of other materials such as microporous 
systems |12| and proteins in hydrated state |13) . The reported values of the coefficient p vary significantly, 
about three orders of magnitude from less than one up a few hundred [U [2 [3j IH [5] . However, for a variety 
of materials, the p value falls mostly in the range from 0.5 to 10 [13], while in the majority of the cases the 
p value is near unity, as it is shown as well in Fig. 3 of Ref. [S]. 

The factors affecting the accurate estimation of p value are discussed in Refs. [HI [TS]. However, some 
aspects should be mentioned here which are related to the electrode effects. In the case of fully-blocking 
electrodes, if the data do not extend to the region where cr'(a;) decreases towards zero in the lower-frequency 
plateau, the result might not define the DC conductivity accurately. In the case of partial-blocking electrodes 
behavior two regions of constant (j'{uj) values could possibly appear [TB]. In such cases it could be possible 
that the higher-frequency plateau region would lead to a more plausible p estimate than the lower-frequency 
one, even though the latter is considered as the DC conductivity [16) . 

In the various models, which have been proposed for the description of the dielectric response of disordered 
conductive materials, the BNN relation has been used as a testing equation through the calculation of p 
coefficient [71 El [171 [18]. The value of p is definitive, in order one to classify conductive materials. According to 
Hunt [19], the role of Coulomb interactions in the derivation of the BNN relation is of great importance, while 
the non-universality of the high frequency limit of the AC conductivity is incompatible with a universality 
in the BNN relation [20]. However, various models proposed for a particular kind of materials such as ionic 
glasses and disordered conductors, give universal values for the BNN coefficient p. Dyre [7] obtains a value of 
p = 0.42 by using the random energy barrier model. In a subsequent work Dyre and Schr0der [S] reported a 
value of p = 1.5±0.4 for the simulation of the symmetric hopping model in the extreme disorder limit. It has 
been pointed out by Macdonald [5] that the Kl conducting-system model could lead to a quantitative value 
for p, which depends on the value of /3ic of the Kohlrauch-Williams-Watts (KWW) stretched exponential 
response function. For ion-conducting homogeneous glasses and single crystals with charge motion allowed 
in all three dimensions, it has been shown that the only possible value is f3ic = 1/3 and the resulting high 
frequency limiting response power law exponent is equal to 2/3 [18]. According to these values, the BNN 
coefficient has a universal value oi p = 1.65, while in the framework of the Kl model, as /3ic — >■ 1, p should 
also approach unity in the limit. In a recent paper Macdonald [T5] has provided a detailed analysis and p 
estimates for the variety of conductive-system models. These models involving a single fractional exponent, 
for an appreciable range of exponent values, show that the p values are quite near 1. 
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A modified BNN relation has been suggested by Dygas in [21] . It is proposed that the values of modified 
P coefficient are related to the spatial extent and time scale of nonrandom local hopping of charge carriers. It 
also gave an expression of the BNN coefficient in the case of Cole-Cole dielectric behavior of ionic polarization 
mechanism. 



The BNN equation quantifies the relation between short range and long range ions motion of the AC 
response of conductive materials. To be specific the p coefficient reflects a measure of the correlation between 
AC and DC conductivity. In the present work, we will attempt to derive the BNN coefHcient based on 
impedance spectroscopy formalism, as well as on widely used phenomenological and empirical relations and 
to discuss the results with relevant published works. The exact knowledge of the parameters on which the 
p coefficient depends, is of great importance, not only from the fundamental point of view, but also for 
applications, because this could lead to design and development of a variety of materials with predetermined 
dielectric and electrical properties. 



2 Theoretical considerations 



The complex conductivity cr*(w) — cr'^u) + ja"{uj) is connected to the total complex dielectric constant 
e*{uj) — s'{uj) — je"{uj) via the following relation 

a*{Lj)=jujeoe*{uj). (2) 

In the above relation, if the contribution of the DC conductivity, ctq, is subtracted from s*{uj), then 

a*{uj) = (JQ + jujeo£*^{uj) , (3) 

where s'^iuj) represents the complex dielectric constant caused only from the dynamic conductivity. 



In ionic materials, the description of the real part of complex conductivity spectra in the low frequency 
regime, below 100 MHz, and in the absence of electrode polarization effects, is given by the equation [22l [23l 



a'{uj) = CTo 



l+(- 



(4) 



where n (with < n < 1) is a constant. The characteristic frequency ujq corresponds to the onset of AC 
conductivity. At this frequency, the real part of complex conductivity becomes twice to that of the DC 
conductivity, a'{u!o) — 2ao. This last equation has been introduced in Ref. [27], in order to describe crystals 
with defects and an activated number of charge carriers. Eq. (|4|) cannot be taken as a model relation, because 
it is not able to reproduce the individual AC response characteristics which are derived from other functions of 
impedance spectroscopy (i.e. the peak and the dielectric strength, in s*(ui) formalism when ionic dispersions 
take place). Eq. is considered as a relation which approximates well only the frequency dependence of 
the real part a' , since at low frequencies, a' describes the DC conductivity plateau of AC response, while at 
high frequencies, a' describes the well-known Jonscher power law behavior [28]. In general, depending on 
the individual characteristics of s'^{uj) responses, Eq. (j4|) can or cannot describe satisfactorily the function 
cr'(w) at the onset region. 
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If we assume that the real parts of Eqs. ©, ^ and Eq. ^ are equal not only at the high frequency 
hmit but also at uj = ujq, then in these cases, from Eqs. ([2]), ([3]) and (|4]) with cj = cjq wc find 

(To = £otJoe"(wo)/2 (5) 

and 

(To = so^^Qs'di^a) ■ (6) 

So, in these cases the characteristic frequency uiq, should be defined also as the frequency at which the 
losses from the DC conductivity are equal to the respective ones of the dynamic conductivity, since e"(aJo) = 



3 The BNN relation 



In what follows, let us consider that, the conductive response of a disordered material in the frequency 
spectrum under study, is characterized only by the existence of DC conductivity and an AC conductivity term. 
The latter is considered that includes entirely contribution due to mobile ions effects. These effects should 
lead to the appearance of a polarization mechanism in e* formalism, which should take place around the 
onset frequency wo, with loss peak frequency, ujmax, and strength Ae. The real part of complex conductivity 
should be given by using Eq. ([3]), as follows 

a' (uj) ^ ao + eouje'^iuj) . (7) 



The s'^{u)) should be considered as s"{uj) conductive-system values, while for its description the well 
known and widely used Havrilak-Negami (H-N) empirical dielectric function is used here as well |29j . 

, Aesm{l3<j)) 
^d\^) - 



[1 + 2(c^/t^ffjv)" cos(a7r/2) + {oj/ujHN)^''f^^ 

(8) 

where 

= ,,etan ( H^hnT Man/2) \ 
^ Vl + (^/^ffw)"cos(a^/2)y ^ ' 



The shape parameters take values in the range < a, /3 < 1 and are closely related to the slopes in log e'J 
vs. logo; plots (e^' ^ at w ^ uJmax and e'J ^ uj~"^ at w ^ ojmax)- It should be mentioned here, that there 
exist models which have been proposed for the interpretation of the limiting behavior of the H-N relaxation 
function, a behavior that is related to the slopes a and —a/3 [501 131] ■ The limiting case oi a = /3 — 1 
corresponds to Debye behavior with a single relaxation time r = 1/ujhn- The frequency ujhn is related to 
UJmax through 

UJmax = AujHN , (10) 
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where 

/ sin(a^/(2/? + 2)) y/" 

^siii(a/37r/(2/3 + 2))y ' ^ ' 



The total dielectric losses should be written as 



'(^) , 



(12) 



where £"(a;) stands for the losses due to the DC conductivity, ctq- So, from the real part of Eq. (jSj we get 
CTo = £o'^£c('^)- frequency u — loq and with the assumption that the real part of Eqs. ([3]) and ^ are 
equal a,t uj — ujq, the relation £"("^0) — ^di^a) holds (see Eq. ^). Therefore the relation o-q = ^o^o^di^o) 
will be valid and with the help of Eqs. (|8]) and (|9]) we arrive at the BNN relation (Tq = peo^S'^maa: where the 
BNN coefficient is equal to 



P 



Wo /sin(a/37r/(2/3 + 2)) 



1/a 



sin(a7r/(2/3 + 2)) 

sin(/3(^) 



(13) 



where in this last expression 
takes its simplest form, — a,Tctaii{uJo / uj h n) ■ 



[1 + 2(wo/wjTjv)" cos(Q7r/2) H 
is given through Eq. ([9|) at w 



- ujQ, while in Cole-Davidson case (a = 1) 



At the high frequency limit, cu » ojhn, Eq. (jS]) is given by 

e'^iuj) - Aesin(a/3^/2)(c^/wffjv)-"'' , 
while Eq. (jH]) becomes 4> = a7r/2 and Eq. ([7]) is written as 

(t'{uj) ^ (To + £oWHArA£sin(a/?7r/2)(w/a;HAr)i~"'^ . (14) 

The function a'{uj) as given by Eqs. (|4|) and ([M]). in the high frequency limit of these two equations, 
u) ^ t^HN, should be the same. This means that the power law exponent of high frequency limit is n ~ 1 — a/3, 
while 

n ^0 n-1 r\ 



" eoA£sin(a^7r/2) 

holds as well. The previous relations Eqs. ([HI) and (|15p . coincide with the ones in Ref. [21], where similar 
expressions have been derived previously. Then based on Eq. ([T5|) and with the use of Eqs. ([T]), ([TU| and 
(ITT]) the BNN coefficient should be expressed as 

wo . (al3Ti\ /sin(a/?7r/(2/3 + 2))^ 



^=l^j ^"^l— jUin(W(2/3 + 2)) j ' 

This last expression, in the Cole-Cole case, coincides with Eq. (15) of [21 . Eq. (IT51) of this present work 
could be extracted as well from Eq. (15) of Dygas work [21] with the help of our Eq. ([TT|). However, it 
seems that an equation like our Eq. ([161) ^lot been reported there, since the research work in Ref. [21] is 
concentrated between other things in the physical meaning of their modified BNN coefHcient. 
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Figure 1: The BNN coefficient p as a function of H~N parameter /3 for different frequency ratios 
(10, 20, ... , 160) according to Eqs ([T3)) . (fT6)) . This figure represents the C-D behavior. 



4 Discussion 



A distinction of the analysis presented above for the extraction of both Eqs. (|13p and (fT6)l . in relation to other 
conductive-system dispersion models which are deahng with BNN matters, is the use of the H-N dielectric 
function here. The H-N function should be considered as e"(a;) conductive-system values. The fact, that 
the H-N function is included in the AC conductivity term, is a consequence of its definition. The coefficient 
of the BNN relation has been found to depend on the H-N parameters a, /3 and the frequency ratio luq/ujhn 
in both Eqs. (fT3|) and ([16]). The parameters a and (3 represent the width and the skewness of the dielectric 
loss, s'^{uj), when viewed in a log(w) plot, and also describe the distribution function of the relaxation times 
[32) . The ratio ujq/ujhn should be viewed equivalently as the ratio of two characteristic time lengths, tuN /to, 
where tuN is the characteristic time length of the ionic polarization mechanism and to is the respective one 
of the onset of AC conductivity. 

In Figs. 1 and 2 the dependence of BNN coefficient p versus the shape parameters a,/3 is given, for 
various values of the frequency ratio loq/lohn according to Eqs. (jl3l) and (jl6p . Two characteristic cases could 
be observed there: the asymmetric Cole-Davidson (C-D) behavior (Fig.l) and the symmetric Cole-Cole 
(C-C) one (Fig. 2). As it is shown in both figures, for a behavior close to Debye one, the BNN coefficient 
approaches 1. For particular values of the parameters a, /3 in both cases, the value of p increases by increasing 
the frequency ratio. 
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Figure 2: The BNN cocfBcient p as a function of H-N parameter a for different frequency ratios 
(10, 20, ... , 160) according to Eqs ([T3)) . (fT6)) . This figure represents the C-C behavior. 



In both cases (C-D and C-C behavior) the value of p increases up to a maximum value as /3 or a decrease, 
and then p decreases as /? or a continues to decrease. Values of p < 1 are allowed in Eqs. (fT3|) and (fT6| for 
both behaviors. We observe also that for the majority of /3 or a parameter values, and for the ratios ljq/luhn 
shown in Figs. [T] and [2] the p coefficient takes values between 0.2 and 10. As the ratio ujq/ojhn gradually 
decreases, the p coefficient tends to take values near unity. 

As it can be seen in Fig. [l] the coincidence of the p coefficient values as specified by Eqs. (IT3)) and is 
excellent. Two different functions have identical behavior in C-D behavior; only for lower frequencies ratio 
ujq/ujhn is there little difference in the values of p, and this difference is less than 0.07. It is also apparent 
in Fig. [2] in the C-C behavior, that the values of p coefficient as provided by Eqs. (|13|) and (fT6|) differ. But 
both p values are similar to the change as well the a parameter changes. The values of p coincide with each 
other to a greater range of high values of a parameter, as the frequency ratio increases in C-C case. 

To clarify the discrepancy in symmetrical C-C behavior, it is necessary to present a representative 
simulation. For parameter values a = 0.4, /3 = 1 and frequency ratio luq/whn = 40, the values of BNN 
coefficient according to Eqs. ([T3l) and (fT6| are found to be 3.78 and 5.38 respectively. Choosing the values 
(To = 10~^ S/m, uiHN = i^max — 100 rad/s and the requirement to satisfy Eq. ([ij, the extracted values are: 
Ae = 299 from Eq. ^ and Ae = 210 from Eq. p^ . 

The corresponding curves are presented in Fig |3] with different Ae values of Eq. (O (with Eqs. ([8]), ([9])). 
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Figure 3: Simulation curves of the AC conductivity a' as a function of w (details in the text). 

Also presented the corresponding Eq. ^ with ctq — 10^^ S/m, wq = 4000 rad/s and n = 0.6 is also presented. 
It is clear that only the result of Eq. (fT6|) produces a curve based on Eq. ©(with Eqs. ^ and ([9])), which 
is described very well by the Eq. (j4]) at low and high frequencies limits. The resulting Ae as calculated 
using Eq. (fT3|) fails to describe the power law frequency dependence of AC conductivity in the case of C-C 
dielectric behavior. 

It should be noted here, that the additional assumption made in order to extract Eq. (IT51) is, that Eq. ^ 
and Eq. ([7]) (with Eqs. ([U, ^) should be equal (or at least approximate to a great extent), not only to 
high frequencies limit, but also at the characteristic frequency wq- The failure of Eq. (fT3| in C-C dielectric 
behavior could possibly due to the fact that the previous additional assumption cannot be applied in these 
symmetric responses. As it is obvious in Fig. [3] the power law frequency dependence of AC conductivity as 
given by Eq. ^ differs substantially at a; = wo, with the one given by Eq. ([7]) with the aid of Eq. (|16p . 
In contrast, Eq. ([T]), based on the parameters of Eq. (jlSp . although is not approximated well with power 
law (Eq. (|4])) at the high frequencies limit seems to be approximated at frequency lo = ujq very well. The 
identification of p values in the asymmetric C~D behavior suggests that in these cases, the function <j'{uj) 
approximated well by Eq. Q not only at the high frequency limit but also at frequency lo — loq. This should 
be the reason that the relations Eq. and Eq. although different, actually coincide in the prediction 
of BNN coefficient, in C-D behavior. 

In what follows trying to achieve contact with real materials, we will refer to a case study of dielectric 
response of a crosslinked polyurethane (PUR) which satisfies the BNN relation. The e"{uj) experimental 
data of PUR ii^ are fitted with a sum of H-N expression (Eqs. ([8]) and ([9])) and the term Alo^^ which 
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represents the contribution of the DC conductivity, ao, to the dielectric losses spectra with value of fc = 1. A 
representative fit is shown at a temperature 363 K in Fig. |4l The dielectric dispersion exhibits, not only at 
this temperature, C-D behavior. The values of ctq where resolved from the extrapolation of low frequencies 



£ 




10 



,-1 







f(Hz) 



Figure 4: The imaginary part e" of the complex dielectric constant as a function of of frequency / at a 
temperature 363 K for a crosslinked polyurethane. In the figure is shown the total fitting as well as the 
contribution of the polarization mechanism and the DC conductivity according to fitting process are shown. 

(T'(a;) plateau at 1 mHz. The /o values (wo — 27r/o) were taking these corresponding to values 2aa according 
to a'^uj) data. It should be noted that the dipolar dispersions appear at higher enough frequencies [33 and 
so their contribution to s"{uj) and a'{uj) frequency spectrum studied here is negligible. Also during fitting 
process of e"{uj) data, as well as the estimation of the DC conductivity, ctqi we have been restricted at a 
suitable low frequencies region, where the influence of the electrode effects could be considered as negligible. 
The parameter values resulting from the fitting process in e"{Lo) data, as well as the values of co ^-nd /o 
which have been calculated from a'{uj) data are given in Tabled] 

Experimental data of ct'(w) are presented in Fig. [5] In the same figure, the simulation curves of Eq. (|4]) 
and Eq. (O (with Eqs. ([8|), (O) are shown according to the parameters values of Table [TJ Despite the fact 
that Eq. dU does not describe adequately the behavior of the rest two curves near the onset region of AC 
conductivity, there is a very good match in all curves for frequencies uj > loq ■ 

The values of the coefficient p, as calculated based on the predictions of Eq. (IT^ and Eq. ([TC]) as well as on 
the BNN relation, Eq. ([T|), are presented in Table[2] In order to do these calculations the corresponding values 
contained in Table [T] have been used. Although the temperature window is relatively small , it is evident 
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Table 1: The parameter values of a, /3, Ae and /hn resulting from the fitting process in s"{uj) data as well 
as the values of (Tq and /o for crosslinked polyurethane at different temperatures T. We consider that the 
values of ctq and /o have been estimated with high accuracy and therefore those are given without errors. 



T(K) 


a 


/3 


Ae 


fHN (Hz) 




/o (Hz) 


o-Q (S / m 


) 


323 


1 ±0.00 


0.39 ±0.01 


4.9 ±0.1 


62 ±5 


132 ± 11 


560 


3.50 X 10" 


-8 


333 


1 ±0.00 


0.37 ±0.01 


4.8 ±0.1 


123 ±5 


272 ± 11 


1185 


7.00 X 10- 


-8 


343 


1 ±0.00 


0.34 ±0.01 


4.9 ±0.1 


198 ±7 


470 ± 17 


2178 


1.27 X 10" 


-7 


353 


1 ±0.00 


0.33 ±0.01 


4.8 ±0.1 


367 ± 11 


894 ± 27 


3560 


2.15 X 10" 


-7 


363 


1 ±0.00 


0.33 ±0.01 


4.7 ±0.1 


656 ± 20 


1597 ±48 


5501 


3.48 X 10- 


-7 



that p coefficient shows a gradual reduction, with increasing temperature, in all calculations included in 
TableO The change of p, throughout the temperature range, appears to be 0.15-0.19. At each temperature 
the p values as they have been calculated using the three equations, show an excellent agreement at all 
temperatures According to this analysis, in order for the coefficient p to be temperature independent, the 

Table 2: The values of p coefficient at various temperatures. 



r(K) 


p (Eq. m) 


p (Eq. m) 


P (Eq. ©) 


323 


1.04 ±0.05 


0.97 ±0.06 


0.98 ±0.07 


333 


1.03 ±0.03 


0.97 ±0.03 


0.97 ±0.04 


343 


1.05 ±0.03 


0.99 ±0.03 


0.99 ±0.03 


353 


0.93 ±0.02 


0.88 ±0.02 


0.90 ±0.03 


363 


0.85 ±0.02 


0.79 ±0.02 


0.83 ±0.03 



parameters a, /3 and ujq/lohn must be temperature independent, or in the case that they are temperature 
dependent, their contributions on p value must be such that they cancel each other. In any other case the 
coefficient of p is expected to be temperature-dependent. 

In Ref. |21| . experimental data have been mentioned, which are relevant to Havriliak-Negami dielectric 
function and the characteristic frequency ratio in three different systems which satisfy the BNN relation. The 
BNN coefficient is calculated according to original BNN relation, and it has been found that it is equal to 
7.1, 1.34 and 2.1 for a single crystal, a type of glass and a PEO polymer electrolyte respectively (for details 
see Table 1 of Ref. [2T). Based on the values of Table 1 of Ref. [21], the BNN coefficient was found to be 
4.6, 1.15 and 1.7 by using Eq. (IT51) . while by using Eq. it was found to be 4.8, 1.32 and 2.1 for a single 
crystal, glass and PEO polymer electrolyte respectively. We observe that the values found using Eq. are 
identical to those based on the original BNN relation for glass and PEO polymer electrolyte, while those of 
Eq. result in values close enough as well. However, in the case of the single crystal, although Eq. ([T5)) 
and Eq. (1161) result in very close values, these values differ significantly from the one extracted using the 
original BNN relation. This difference is possibly due to uncertainties during the fitting process. The fitting 
of Havriliak-Negami function in e"(w) and e'{uj) data after subtraction of the DC conductivity and electrode 
effects respectively, lead generally in scattering of the data points at frequencies lower than the ones which 
correspond to the peak of the dielectric dispersions. This fitting process has as a result uncertainties in the 
estimation of the Havriliak-Negami parameters. 
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Figure 5: The AC conductivity a' as a function of frequency / at 363 K for a crosslinked polyurethane. The 
hues correspond to simulation curves according Eq. (|4]) and Eqs. (O, (|8]), (|9]), and using the corresponding 
values of Table [TJ 

As can be seen observing the values of Table [U the relation in Eq. ([T5|l is satisfied with sufficient accuracy. 
Note here that the strength of this relationship requires the power law Eq. (U) to approach very well the 
function g' (uS) only at the higher frequencies limit, and thus its validity is more general . From Eq. (fT3|) and 
Eq. ([T6| is evident that there are many possible -p values. Also, for a given value of BNN coefficient, there 
are several possible combinations of parameters a, /3, and ujf)/LJHNi according to Eqs. (|13p and ()16p . In these 
cases, not only the high frequencies slope of function (7'(aj), but also the ratio o'o/'^o^e could be different 
with consequences in scaling diagrams of AC conductivity. Of course one question is if all these different 
combinations of parameters may correspond to the response of the real materials. Assume for simplicity 
that the ionic polarization mechanisms of a category of materials, characterized by particular values a, /3 
and loq/ujhn- Then, these materials characterized by a single value of p. This seems to be compatible with 
the universal value of p, as proposed by Macdonald [18^ for ion-conducting glasses and single crystals. His 
conclusion is based on the observation that these materials are characterized by the same shape parameter 
Pic of KWW response function and also the same characteristic frequencies ratio of Kl model. This indicates 
that H-N and KWW shape parameters are closely related, as well as the characteristic frequencies ratio of 
our approach {ujq/ujhn) to that of the Kl model. Particular values of the parameters a, (3 and loq/ujhn 
implies that the corresponding materials not only have same p value but also have the same ratio CTo/i^oAe, 
according to Eq. ()15p . Of course there are many possible combinations of ctq , i^o and Ae values which keep 
this ratio constant. But only some of these combinations correspond to the response of real materials. 

In a recent paper Dyre et al. | 22| applied the fluctuation-dissipation theorem in ionic conductors. 
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According to their calculations, the DC conductivity is 



where q denotes the ions charge, is the Boltzmann constant, T is the absolute temperature, ra(io) is the 
number density of mobile ions at time = I/ujq (loq being the onset frequency), (Ar^(to)} is the single- 
particle mean square displacement, while 7 and H are numbers which are roughly of order one. Roling et 
al. |34j using linear response theory have found a relation between the ionic dielectric strength, Ae, and the 
mean-square displacement of the center of the mobile ions, (i?^(oo)), in ionic glasses conductors 

eoAe = ^{R\oo)) , (18) 

where Ny is referred as the number density of mobile ions. Similar expressions to that of Eq. (jlSp have 
already been reported in the literature [131 |3S] . Now, substituting the quantity ojq / ctq as given by Eq. ([T7)) 
and that of EqAs as given in Eq. ([T5|. back into Eq. (|15p . we arrive at the following expression 



^0 ^HN 



^HN to 



7iJsm(— -)- 



n{to) (Ar2(to)) 



(19) 



The factor 7iJ sin(a/37r/2) in the previous equation takes values close to unity, for typical broad ionic dis- 
persions. The characteristic frequencies or times ratio depends on the number densities of mobile ions at 
different time scales which characterize the AC response, as well as to the corresponding distance lengths 
of ions motion. The other factor affecting the characteristic ratio is the product af3 which represents the 
high frequency absolute slope of ionic dispersion in loge^' vs. logw plot, and it is connected to the high 
frequency limiting response power law exponent of a' . The lower the value of a/3, the stronger its influence 
on the characteristic ratio is. As a consequence of this, different values of the power law exponent in the 
high frequency limit of AC conductivity, lead to different values of BNN coefficient, in accordance with the 
conclusions about the non-universality of p value, as has been suggested by Hunt .20 . From Eq. ([T9l) . which 
is based on two different models [22l|34], it is obvious that particular values of a, (3 and ojq/ujhn correspond 
to a specific correlation of the number densities as well as of the distance length of mobile ions at different 
time scales, which affect the macroscopic characteristic of long and short range ions motion. 



Finally, it should be noted that the high frequency slope of cr', is n = 1 — f3KwW: where Pkww is 
the corresponding KWW stretched exponential parameter [35JEZ|. So, according to the present analysis 
Pkww ~ where a and (3 represent the H-N shape parameters of the ionic polarization mechanism in e* 
formalism. The connection between the H-N and KWW parameters appears also in other relevant works in 
the literature [5^1551 155]. 



5 Conclusions 



Summarizing this work, the BNN relation has been extracted while it has been found that p coefficient 
depends on the H-N shape parameters a and (3 as well as on the ratio ujq/ojhn of two frequencies which 
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are directly connected to characteristic time scales of ions motion (Eq. p^ ). This equation is compared 
to the BNN relation (Eq. (fT6|) ) which could be extracted based on Dygas findings in his previous work. 
Both approaches highlight the dependence of BNN coefficient from the same parameters. In C-D dielectric 
behavior of ionic polarization mechanism, both relationships are nearly identical predictions with very high 
accuracy for the p coefficient, while deviations are observed in the C-C one. Representative simulations 
have shown that only Eq. (jl6p . is consistent with the power law frequency dependence of AC conductivity in 
symmetrical C-C cases. The identification of p values in asymmetric C-D behavior indicates that in these 
cases, the function (t'(w) should be approached very well from the power law not only at high frequency limit 
but also at a; = wq. Both expressions of BNN relation could result in many possible p values. However, the 
response of matter could exhibit individual characteristics. In both cases, of C-D and C-C behavior, the p 
values, depending on ujq/ujhn ratio, lie within the range where the majority of the p values recorded in the 
literature falls. Specifically, p values close to unity arise for the lower values of the ratio loq/ujhn and for an 
appreciable range of the H-N shape parameters /3 and a. 

The two expressions of BNN coefficient, applied to experimental data of a crosslinked polyurethane, which 
satisfies the BNN relation and the ionic polarization mechanism exhibits C-D behavior. Both predictions 
are in excellent agreements with the values calculated from original BNN relation, according to experimental 
data analysis. The results also showed a trend of gradual reduction of the p coefficient, as the temperature 
increases. 
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